function meta_mismatch_plots(proj_meta,cur_site,cells)
all_act={};
pl_cells=[];
mean_fluo=[];
for tnd=1:size(proj_meta(cur_site).rd,2)
    cnt=0;
    for vnd=1:size(proj_meta(cur_site).rd,1)
        all_act{tnd}(cnt+1:cnt+size(proj_meta(cur_site).rd(vnd,tnd).act,1),:)=proj_meta(cur_site).rd(vnd,tnd).act;
        mean_fluo(cnt+1:cnt+size(proj_meta(cur_site).rd(vnd,tnd).act,1),tnd)=mean(proj_meta(cur_site).rd(vnd,tnd).act,2);
        cnt=cnt+size(proj_meta(cur_site).rd(vnd,tnd).act,1);
        pl_cells(vnd)=cnt;
    end
end
fb_snip={};
pb_snip={};
fb_av={};
pb_av={};
for tnd=1:size(proj_meta(cur_site).rd,2)
    nbr_frames=proj_meta(cur_site).rd(1,tnd).nbr_frames;
    stim_id=proj_meta(cur_site).rd(1,tnd).ps_id;
    runP=proj_meta(cur_site).rd(1,tnd).velP_smoothed>0.01;
    runM=proj_meta(cur_site).rd(1,tnd).velM_smoothed>0.01;
    runM_raw=proj_meta(cur_site).rd(1,tnd).velM_smoothed;
    if nargin==3
        for lnd=1:length(cells)
            act=all_act{tnd}(cells(lnd),:);
            [fb_av{lnd,tnd},pb_av{lnd,tnd},fb_snip{lnd,tnd},pb_snip{lnd,tnd}]=pert_ana_meta(nbr_frames,stim_id,runM,runP,runM_raw,act);
        end
    else
        for lnd=1:size(all_act{1},1)
            act=all_act{tnd}(lnd,:);
            [fb_av{lnd,tnd},pb_av{lnd,tnd},fb_snip{lnd,tnd},pb_snip{lnd,tnd}]=pert_ana_meta(nbr_frames,stim_id,runM,runP,runM_raw,act);
        end
    end
end

if nargin==3
    for znd=1:length(cells)
        meta_analysis_LFMplot2(fb_snip,pb_snip,fb_av,pb_av,znd)
        pause
    end
else
    for znd=1:size(all_act{1},1)
        meta_analysis_LFMplot2(fb_snip,pb_snip,fb_av,pb_av,znd)
        pause
    end
end